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Abstract 

Through a complex shift of the time coordinate, a modification of Bjorken flow is in- 
troduced which interpolates between a glasma-like stress tensor at forward rapidities and 
Bjorken-like hydrodynamics around mid-rapidity. A Landau-like full-stopping regime is 
found at early times and rapidities not too large. Approximate agreement with BRAHMS 
data on the rapidity distribution of produced particles at top RHIC energies can be achieved 
if the complex shift of the time coordinate is comparable to the inverse of the saturation 
scale. The form of the stress tensor follows essentially from symmetry considerations, and it 
can be expressed in closed form. 
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1 Outline of the construction 



The central idea of Bjorken flow pQ is that there is approximate boost invariance in the 
direction of the beamline for the dynamics near mid-rapidity. Particles that form near mid- 
rapidity are assumed to do so at some definite, boost-invariant proper time r form m lfm/c, 
where r = \/t 2 — x\. Assuming that all motion is in the beamline direction, the individual 
four-velocities of the produced particles can be deduced to have the form 
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To the extent that one may use a hydrodynamical description, the boost invariance also 
constrains the local energy density e: it can depend on t and £3 only through the combination 
t. If one assumes translational and rotational invariance in the collision plane, then e cannot 
depend on x\ or x<i- Thus e = e(r). One can obtain explicit expressions for e(r) by assuming 
a specific form for the stress tensor. For example, let's choose 

T^v = EUfj,Uu + -{g^v + u^Uu) , (2) 

corresponding to inviscid, conformal hydrodynamics. Then one immediately finds 
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where eo is a constant. 

A trivial modification of Bjorken flow is to set 



u? = I , ' "° =,0,0 
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where t3 is a constant. If is real, then we have simply translated Bjorken flow in time. If 
t3 is complex, we have something new, and all the quantities with a superscript C become 
complex. The complexified stress tensor T^ u still obeys the conservation equations V^T^ = 
0, and because these are linear equations, a conserved stress tensor with all components real 
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can be obtained by denning 

T, v = KeT^. (5) 

as identified in ^ will not in general satisfy the inviscid hydrodynamic constitutive 
relations. 

It may seem that @-(§ are an unmotivated and unpromising line of attack on the 
problem of describing the rapidity structure of heavy ion collisions. However, provided the 
phase of e c is chosen correctly (namely arge c = 7r/3 when argt3 = vr/2), an appealing 
spacetime picture emerges in the forward lightcone, as illustrated in figure [Tj At early times 
and rapidities not too large, there is a full-stopping region reminiscent of the Landau model. 
Although the hydrodynamic constitutive relations hold to good accuracy in this region, 
they do not hold uniformly in its causal future. Instead, one recovers the hydrodynamic 
constitutive relations, and Bjorken flow, asymptotically at late proper times r with rapidity 
held fixed; but at forward rapidities one obtains a glasma-like form of the stress tensor, with 
longitudinal pressure nearly equal to minus the energy density. Interestingly, at r = |t 3 |, 
there is a very simple relation between the rapidity yp of the fluid and the pseudorapidity: 

77 

y F = - when r = |t 3 | . (6) 

This is to be compared with the relation yp = r\ for Bjorken flow. 

The organization of the rest of this paper is as follows. In section [2] I explain how the 
phase of e c is fixed by positive energy considerations. In section [3] I explain the main features 
of the rapidity structure of the flow, providing in particular a more quantitative version of the 
diagram in figure [TJ In section [4] I consider a simplified version of hadronization and exhibit 
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a comparison of the predicted rapidity profile with data from the BRAHMS experiment. I 
conclude in section [5] with a discussion of symmetry groups and local entropy production. 



2 Positive energy constraints 

In addition to the conservation equations, a sensible stress tensor must satisfy some sort of 
positive energy constraints. The conditions we will apply are that the Landau frame can be 
defined throughout the future light-wedge of the collision plane, and that the energy density 
in the Landau frame is positive throughout this wedge. 

To pass from the laboratory frame to the Landau frame, we must apply a boost which 
makes the stress tensor diagonal^] This boost is in the x% direction, because the only non- 
diagonal components of T^ u are T 03 = T 30 . Suppressing the X\ and x 2 directions for brevity, 
we may express 

T _ f T oo T 03 \ _ fcoshy F sinh?/ F 

\To3 T 33 J \smhy F coshyp 



Tp V — A%A^T a/3 
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where e L is the Landau frame energy density and p 3 is the Landau frame longitudinal 
pressure. The transverse pressure Pt = Pi = P2 is the same in the Landau frame as it is in 
our original frame. It is easy to show from ^ that 

1 1-1 2T 03 
Uf = — tanh 



€ 



2 Too + T33 

l _ Too — T33 + a/ (Tqo — 2T 3 + T 33 )(Too + 2T 3 + T 33 ) 



l ~ T 00 + T33 + a/ (T o — 2T 03 + T 33 )(T o + 2T 03 + T 33 ) 
P3 = " o 



The Landau frame exists precisely when yp as defined in (|8j) is real. I claim that in order for 
the Landau frame to exist throughout the future light-cone of the collision plane, one must 
set 

$ = * 2/3 $, (9) 

where is real. Furthermore, e L > everywhere in the future light-cone provided > 0. 
A full demonstration of these claims is tedious, but I will give sufficient indications here to 



1 For a more careful discussion of how to define Landau frame, see for example [2]. 
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show that no other phase than the one in ^ will suffice. To this end, let's set 



£o c = z 2 /3 e «> tg = i T= ^J t 2_ x 2 = 1} (io) 

where without loss of generality we can limit arg# G (— 7r/2,7r/2) by allowing e R to be a 
signed real quantity (anticipating that the constraint e R > will emerge from later argu- 
ments). Straightforward computation now leads to 

tanh2?/ F = - — — — = 1 + — — , (11) 

-?oo + T 33 t \ t cos 9 + (t 2 - 1) sin 6/ / 



It should be borne in mind that (11) holds only with the special values (10), in particular 



1. If sin# ^ 0, then by expanding (11) at large t one can see immediately that the 



right hand side is greater than 1 when t is sufficiently large. Thus 9 = 0, and (11) simplifies 
dramatically to yp = r]/2, where 

r] = tanh" 1 — (12) 

is the spacetime pseudorapidity. This result, previously quoted as ([6]), contrasts with the 
result jjf — V for standard Bjorken flow and provides some preliminary indication that the 
flow is more focused near mid-rapidity. It is also useful to note that yp = when rj = for 
all t, and upon setting 9 = one finds 

^2/3 £ R 

Too = e L = 3p^ = 3p T = Re - ° when x 3 = . (13) 

{t + hr 6 

So the ideal hydrodynamic constitutive relations are satisfied precisely at mid-rapidity, and 
^oo ~ e R /2^ 4 ^ 3 & t mid-rapidity for t \ts\, which is the same functional dependence as for 



Bjorken flow. Thus > 0, as promised. 



3 Features of the rapidity dependence 

Up to an overall rescaling of time, t3 = i is the unique choice for the type of flow we 
are interested in. As argued above, up to an overall rescaling of energy, = i 2 ^ 3 is the 
unique choice for a physically sensible stress tensor. Thus, through Q, (pH), and (J8l), I have 
constructed an essentially unique stress tensor. Its main features, exhibited in figures [2](3j 
may be summarized as follows: 

• There is a region at late times in which yp rs r] and pf/e ~ 1/3: Bjorken flow is 
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Figure 2: Top: The fluid rapidity versus spacetime pseudorapidity, at several different values 
of Bjorken time, for t 3 = i. Middle: The Landau frame energy density, scaled by r~ 4//3 , 
versus spacetime pseudorapidity, at several different values of Bjorken time. Bottom: The 
ratio of longitudinal pressure to energy density, in the Landau frame, versus spacetime 
pseudorapidity, at several different values of Bjorken time. The ratio p%/e L is not shown for 
r < 1 because it happens to be identical for r and 1 13 1 2 / t . 
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approximately recovered. To define more precisely where this happens, I note that 



^f>0.3 when U <> -1.4 + log ^ 

e m 

dyp . , 7~ 

~ — > 0.9 when 77 & log , — r . 



(14) 



The energy density and pressure become broader in pseudorapidity, and the energy 
density acquires a characteristic double-hump structure for t <L 10|t 3 |. It is interesting 
that Df attains values slightly larger than r] in the Bjorken region (by roughly 10%). 
Uf then approaches 77 from above as r — ^ 00 at fixed rj. Because of this overshoot, it is 



perhaps better to describe the region where the conditions (14) hold as "Bj or ken-like. 



There is a region at early times in which the fluid rapidity yp is small and the ratio 
p^/e L is close to 1/3. This is reminiscent of full stopping in the Landau model, so I 
will refer to the region in question as the full-stopping region. Because yp = exactly 
at mid-rapidity for symmetry reasons, it is better to use the smallness of dyp/drj (at 
fixed r) to define the full-stopping region. I find that 

^<0.1 and % > 0.3 when - < 0.6 and < 0.2 . (15) 
dr] e L t 3 |t 3 | 

It is probably appropriate to think of the fluid as far from equilibrium in the full- 
stopping region, because although the fluid is close to satisfying the inviscid hydrody- 
namic constitutive relations there, it ceases to do so in much of the causal future of 
the full-stopping region: in particular, in the glasma-like region to be described next. 

At extreme forward rapidities (and for proper times not too small) the fluid rapidity is 
close to satisfying the curious relation dyp/drj = 1/2 — a relation which, as explained 



around (11), is exactly satisfied for all rapidities when r = |t 3 |. In a similar region, I 
find p 3 < 0, with p% f» — e L as the pseudorapidity becomes large. This is the same as 
the stress tensor of purely longitudinal electric and/or magnetic fields, as considered 
in glasma accounts of pre-thermalization dynamics: see for example [3]. At late times, 



^-f < when 77 > 0.2 + log — r 

e ' ' 1*3 1 

< 0.6 when \n\ ^ 0.9 + log ^~ 

drj |t 3 | 



(16) 



6 




Figure 3: Contours of constant p\je L (green) and dyp/drj (red) demarcate regions associ- 
ated with Bjorken-like flow (blue), Landau-like full-stopping (grey) and glasma-like behavior 
(gold). 

A visual summary of the bulleted points above can be found in figure [3j which may be seen 
as a more quantitative version of figure [TJ 

4 A simple version of hadronization 

The presence of a Bjorken plateau in the energy density which widens as one proceeds 
toward late times is, at least qualitatively, a phenomenologically attractive feature. To be 
more quantitative, one should ask what sort of spectrum of hadrons one gets from isothermal 
freezeout. Using a conformal equation of state all the way to freezeout seems excessively 
naive, so let's use the equation of state sometimes referred to as EOS Q [4J. In this model, 

P =( —^- f ^ T>T C , (17) 
where B is the bag constant and T c is the deconfinement temperature. Using the first law 



of thermodynamics one can deduce from (17) that 



7T 2 7T 2 



o,T 4 + B and p = — o,T 4 - B for T > T c . [IS) 
3(T v 9(T 

Here g„ is the number of degrees of freedom in the deconfined phase. EOS Q calls for a 
strongly first order transition at T = T c , and for T < T c one assumes p = 0.15e, which is 
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supposed to describe a Hagedorn spectrum of non-interacting hadronic resonances. A fairly 
realistic choice of parameters is 



B = 0.35 — — ~- #* = 40 T c = 164MeV. (19) 
fm 

It is hard to see how to incorporate an equation of state with a phase transition into the 
stress tensor construction So let's hadronize at the temperature T c , just before the phase 
transition occurs. Usually Cooper-Frye is applied using a lower freezeout temperature, for 
example T^ in m HOMeVj^] The justification for the Cooper-Frye algorithm relies on local 
equilibration prior to freezeout, so it is dubious to apply it outside the Bjorken plateau. I 
will use a freezeout surface running from 77 = — log t^t to 77 = + log t^t . Depending on one's 
precise definitions, this is roughly the extent of the Bjorken plateau. 



By restricting attention to the high-temperature regime (18), we can view the stress 
tensor of EOS Q as that of a CFT plus a positive cosmological constant. It is perhaps 
unsurprising that fluid rapidity and temperature (suitably defined) are the same as for a 
CFT. The remainder of this paragraph is devoted to a more careful demonstration of this 
claim. First note that ordinary Bjorken flow with EOS Q (and no viscosity) has 

e Q = B +T^ = B + £ cft 

eo/3 (20) 
Pq = -B+ -^g = -B + pcFT ■ 



Because in Bjorken flow is completely independent of the equation of state, (20) implies 

(T^ v )q = (T^u)cft ~ Bg^ v . (21) 
Passing to a complexified stress tensor, as in Q, with 

c e ° 

e = B + 771 1 X ^2 2Y2/3 ' ( 22 ) 

((* + ts) 2 - ^3) 7 

one sees that (T^ u )q = (T^) CFT — Bg^, and it follows immediately that the final stress tensor 



ReT^, obeys the same relation: that is, (|21|) applies unaltered. Because the term — Bg 



is frame-independent, the boost required to pass to Landau frame is the same whether this 
term is present or absent. So yp is indeed the same for EOS Q as for a CFT. Temperature is 

2 Hadronization schemes with separate temperatures for chemical and kinetic freezeout typically set the 
chemical freezeout temperature close to 160 MeV, close to our T c . It would be interesting to see whether these 
more sophisticated hadronization schemes lead to significantly different results for the rapidity distribution. 
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a little trickier because the final stress tensor does not satisfy the hydro dynamic constitutive 
relations. The obvious approach is to define 

T = ^bV>\ m 

V.^ 9* J 

This definition follows from plugging the Landau frame energy density into the first relation 



in (18) and solving for the temperature. T as defined in (23) evolves identically to the 
temperature of a CFT, both in Bjorken flow and in complex deformations of it. 

I am going to use a simplified version of the Cooper-Frye expressions for the distribution 
of outgoing particles d 3 N/dp 3 with energy E: namely 

J 3 N g f 1 



E—fY = t/\3 / ^ P Ax^dx u dx p fe u ^ / T , (24) 
dp 6 (27r) d ,/ s 3! 

where = (— coshy^, 0, 0, sinhy^) in the laboratory reference frame and temperature T is 



given by (23). (My convention is to sum over indices without regard to ordering, which is the 
reason we require the explicit 1/3! in the integrand.) The integration surface S is determined 
by the equation T = T freezeout , where T freezeout = T c is the freezeout temperature. The overall 
factor g refers to the degeneracy of produced particles. The rapidity of a produced particle 
can be defined (in the laboratory frame) in terms of its momentum p M as 

y = tanh- 1 1 . (25) 

Standard manipulations lead to an expression for the rapidity distribution of produced par- 
ticles: 



dN f . , „d 3 N 
dy 



where 



and 



f d 3 N f 

J dpidp 2 E -j-p = J (Qrdr + q^drj) (26) 

q T = Q sinh(y - rj) q v = -Qt cosh(?/ - 77) (27) 

gVo\± f°° 2 0Vol_L ^ 2 , 1 2rf\ -ffi , OQ x 

Q = In — vT / « m -L m ± e T = o 9 T + mT + o m T e T • ( 28 ) 

(27r) z y m - 2tt^ V 2 / 



Here Volj_ is the volume in the X1-X2 directions, m± = \jp\ + p\ + m 2 , and 

T = Tsech(y-y F ). (29) 
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(dN/dy) sca i ed 




Figure 4: The rapidity distribution of produced particles in rapidity, normalized to unity 
at mid-rapidity. The grey dots are PHOBOS data for dN/drj, from the lower-right panel of 
Figure 18 of |6]. The black dots are BRAHMS data for dN/dy of positively charged pions, 
from Figure 4 of [7j. The black curve is from the purely hydrodynamic Landau model, as 
explained in [7]. 



In figure [4] I show predictions of dN/ dy as compared to data for central RHIC collisions 
at v /sjvjv = 200 GeV. The main points to note about my parameter choices are: 



1. I consider only pions as the outgoing particles, so I use m = w 140 Me V in (28). 



2. For | -£3 1 , I consider multiples of the length 0.07fm, which is about half the thickness of 
the incoming nucleus. 

3. The energy density at r = 0.7fm and r\ = is set equal to 7.7GeV/fm 3 in order to 
approximate the conditions of a central RHIC collision at y/s^N — 200 GeV (see for 
example [S]). 

4. I restricted the rapidity range of the freezeout surface to \r]\ < logj^y, corresponding 
approximately to the extent of the Bjorken plateau. 

5. The freezeout surface turns out to be fairly close to isochronous (in Bjorken time r), 
with Tfreezcout ~ 2.3 fm. Early freezeout is expected because of the higher-than-typical 
value of freezeout temperature. 

6. I normalize dN/dy to 1 at mid-rapidity. Estimates of the effective Volj_ and inclusion of 
more hadron species would be needed in order to obtain realistic dN/ dy at mid-rapidity. 
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Comparison with the BRAHMS data on dN/dy favors | -t 3 1 pa 0.2 fm, corresponding to an 
energy scale of 1 GeV. It is gratifying that this is close to the saturation scale Q s ~ 1.5 GeV 
for central gold-gold collisions at ^snn = 200 GeV. 

5 Discussion 

The key principle underlying Bjorken flow is boost-invariance in the longitudinal direction. 
Complexification weakens this principle but does not destroy it entirely. To see how sym- 
metries behave in the complexification, let's introduce notation for the generators of the 
Poincare group ISO(3, 1): 

• Translations, T( M ) = d^. Thus, for example, is translation in the x\ direction. 

• Spatial rotations, R(ij) = %idj — Xjdi. 

• Boosts, -B(j) = tdi + Xid t . 

Usually one demands invariance under -8(3). There is essentially one combination of t and 
£3 that is .B( 3 )-invariant, namely r 2 = t 2 — x\. If we instead demand invariance under 

b = B {3) + t 3 T (3) , (30) 

then the invariant combination of t and £3 is (t + t 3 ) 2 — x\. It is straightforward to show that 
JO-brf — 0, where Cb is the Lie derivative and is defined as in Q. Likewise, C^T^ = 0. 
When 13 is imaginary — the case of interest — b is not a member of the algebra of ISO(3, 1), 
but it is in the complexification of this algebra. So the complexified stress tensor is invariant 
under the complexified symmetry. However, the final form ReT^ is not invariant under b 
or any obvious modification of it. 

Although I have focused exclusively on the case c 2 = 1/3, it is worth noting that other 
constant speeds of sound can be treated in essentially the same way. The relation yp = 77/2 
is recovered at t = \t^\ for arbitrary c 2 . Another way to say this is that the construction goes 
through for conformal field theories in dimensions other than 3 + 1. 



Following the definition (23) of temperature in terms of energy density, and restricting 



attention to a conformal equation of state, one may define an entropy current 

s, = T 3 u,. (31) 
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(There should be an overall constant factor on the right hand side of (31 ), proportional to the 
number of degrees of freedom, but this factor doesn't affect the discussion to follow.) When 
the inviscid hydrodynamic constitutive relations hold, the equations of motion V M T^ = 
imply conservation of entropy, V^s^ = 0. The second law of thermodynamics, applied locally, 
requires locally increasing entropy, V^s M > 0. By direct calculation starting from (|4])-([5]), I 
found that entropy locally increases in the Bjorken-like region, but not in the full-stopping 
region, and in only a part of the glasma-like region which is not too close to the lightconej^] 
Failure of the second law outside the Bjorken-like region seems at first alarming; but what 
one should conclude is that (31) is a poor approximation to the entropy current except in 
the Bjorken-like region. This makes sense because only in that region is there good reason 
to think that local equilibration has occurred. 
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